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METHOD AND SYSTEM OF QUANTITATIVELY ASSESSING DIFFUSION 
CHARACTERISTICS IN A TISSUE 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] The present invention claims the benefit of U.S. Ser. No. 60/471,589 filed 
May 19, 2003. 

BACKGROUND OF THE INVENTION 

[0002] The present invention relates generally to a method and system of 
processing diffusion weighted magnetic resonance (MR) data and, more particularly, to a 
tensor-free method and system of quantitative assessment of diffusion characteristics in 
a targeted tissue. 

[0003] When a substance such as human tissue is subjected to a uniform 
magnetic field (polarizing field B 0 ), the individual magnetic moments of the spins in the 
tissue attempt to align with this polarizing field, but precess about it in random order at 
their characteristic Larmor frequency. If the substance, or tissue, is subjected to a 
magnetic field (excitation field B0 which is in the x-y plane and which is near the Larmor 
frequency, the net aligned moment, or "longitudinal magnetization", M z , may be rotated, 
or "tipped", into the x-y plane to produce a net transverse magnetic moment M t . A signal 
is emitted by the excited spins after the excitation signal B 1 is terminated and this signal 
may be received and processed to form an image. 

[0004] When utilizing these signals to produce images, magnetic field gradients 
(G Xl G y , and G z ) are employed. Typically, the region to be imaged is scanned by a 
sequence of measurement cycles in which these gradients vary according to the 
particular localization method being used. The resulting set of received NMR signals are 
digitized and processed to reconstruct the image using one of many well known 
reconstruction techniques. 
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[0005] It is generally well-known that MR imaging techniques may be exploited 
to non-invasively discriminate between normal and pathological tissues in a subject, 
such as a medical patient. One such MR imaging technique is predicated on the 
acquisition of diffusion-weighted MR data that enables a physician, radiologist, or other 
health care provider to assess diffusion anisotropy in a particular tissue region of a 
patient, such as in the human brain. Until recently, health care providers were unable to 
quantitatively assess diffusion. That is, MR imaging supported the acquisition and 
reconstruction of MR data that allowed a clinician to visually inspect and determine, from 
the reconstructed images, a qualitative measure of diffusion. Since assessment of 
diffusion was limited to a qualitative evaluation, the extent of the diagnosis process was 
limited to the experience and expertise of the health care provider as well as the quality 
of the data acquisition process. Moreover, patient movement, field inhomogeneities, and 
the like could negatively affect the data acquisition process thereby hindering a 
physician's or radiologists' ability to accurately assess diffusion when visually 
determining and making a diagnosis. 

[0006] As a result, recently, quantitative approaches to diffusion weighted 
imaging have been proposed. With these quantitative approaches, a value of diffusion is 
determined that enables a clinician to quantitatively assess the diffusion in the tissue 
imaged to determine whether the tissue is normal or whether a pathology is present. 
Moreover, by quantitatively measuring diffusion, the degree of pathology may also be 
determined. 

[0007] However, with some of these known quantitative approaches, a tensor- 
based approach is implemented. These conventional tensor-dependent approaches, 
however, significantly limit assessment of diffusion when diffusion extends or tracts 
along multiple directions in a particular imaging volume. Further, these conventional 
tensor-based techniques must make numerous assumptions regarding the diffusion 
data. As a result, the accuracy of the quantitative assessment is remarkably dependent 
on the assumptions made. 
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[0008] It would therefore be desirable to design a method and system of 
quantitatively assessing diffusion weighted data in diagnosing a tissue without making 
assumptions as required by a tensor-based model that is also applicable to voxels 
having multiple diffusion tracts. 
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BRIEF DESCRIPTION OF THE INVENTION 

[0009] The present invention is directed to a method and system of assessing 
diffusion characteristics in a region-of-interest that overcomes the aforementioned 
drawbacks. 

[0010] The present invention includes a process and system to implement such 
a process that is configured to determine diffusion properties within a region-of-interest. 
The invention may be implemented to determine shared diffusion anisotropy between 
two voxels. Based on the value of a diffusion index or metric, the present invention 
provides an efficient technique of determining if diffusion in neighboring voxels is 
anisotropic and if the diffusion is oriented in the same direction within each voxel. The 
present invention may be used to rapidly determine fractional anisotropy in a single 
voxel or used to filter diffusion weighted imaging data. The present invention may also 
be used to color-code diffusion data such that diffusion and tract orientation are readily 
identifiable in a reconstructed image. 

[0011] Therefore, in accordance with one aspect, the present invention includes 
a method of assessing diffusion of a substance in a region-of-interest. The method 
includes the steps of determining a diffusion value for one or more individual directions 
in at least one voxel and determining a value of average diffusion in the at least one 
voxel. The method further includes the step of, from the individual diffusion values and 
the average diffusion value, determining a diffusion index indicative of diffusion in the at 
least one voxel. 

[0012] According to another aspect, the present invention includes an MRI 
apparatus comprising a magnetic resonance imaging (MRI) system having a plurality of 
gradient coils positioned about a bore of a magnet to impress a polarizing magnetic field 
and an RF transceiver system and an RF switch controlled by a pulse module to transmit 
RF signals to an RF coil assembly to acquire MR images. The system also includes a 
computer programmed to access diffusion-weighted imaging data from a region-of- 
interest and isolate at least two adjacent voxels within the region-of-interest. The 
computer is further programmed to determine a metric indicative of diffusion along 
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multiple directions in the at least two adjacent voxels without fitting diffusion-weighted 
imaging data from the at least two adjacent voxels to a tensor model 

[0013] In accordance with another aspect of the present invention, a computer 
readable storage medium having stored thereon a computer program representing a set 
of instructions is provided. The set of instructions, when executed by a computer, 
causes the computer to determine an individual diffusion value for each of N diffusion 
directions in at least one voxel and determine an average diffusion value for the at least 
one voxel. The set of instructions further causes the computer to determine a diffusion 
index indicative of diffusion anisotropy and diffusion direction in the at least one voxel. 

[0014] Various other features, objects and advantages of the present invention 
will be made apparent from the following detailed description and the drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] The drawings illustrate one preferred embodiment presently 
contemplated for carrying out the invention. 
[0016] In the drawings: 

[0017] Fig. 1 is a schematic block diagram of an MR imaging system for use with 
the present invention. 

[0018] Fig. 2 is a flow chart setting forth the steps of a technique of quantitatively 
assessing diffusion in a region-of-interest. 

[0019] Fig. 3 is a flow chart setting forth the steps of a technique of determining 
fractional anisotropy in a selected voxel. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0020] Referring to Fig. 1, the major components of a preferred magnetic 
resonance imaging (MRI) system 10 incorporating the present invention are shown. The 
operation of the system is controlled from an operator console 12 which includes a 
keyboard or other input device 13, a control panel 14, and a display screen 16. The 
console 12 communicates through a link 18 with a separate computer system 20 that 
enables an operator to control the production and display of images on the display 
screen 16. The computer system 20 includes a number of modules which communicate 
with each other through a backplane 20a. These include an image processor module 
22, a CPU module 24 and a memory module 26, known in the art as a frame buffer for 
storing image data arrays. The computer system 20 is linked to disk storage 28 and 
tape drive 30 for storage of image data and programs, and communicates with a 
separate system control 32 through a high speed serial link 34. The input device 13 can 
include a mouse, joystick, keyboard, track ball, touch activated screen, light wand, voice 
control, or any similar or equivalent input device, and may be used for interactive 
geometry prescription. 

[0021] The system control 32 includes a set of modules connected together by a 
backplane 32a. These include a CPU module 36 and a pulse generator module 38 
which connects to the operator console 12 through a serial link 40. It is through link 40 
that the system control 32 receives commands from the operator to indicate the scan 
sequence that is to be performed. The pulse generator module 38 operates the system 
components to carry out the desired scan sequence and produces data which indicates 
the timing, strength and shape of the RF pulses produced, and the timing and length of 
the data acquisition window. The pulse generator module 38 connects to a set of 
gradient amplifiers 42, to indicate the timing and shape of the gradient pulses that are 
produced during the scan. The pulse generator module 38 can also receive patient data 
from a physiological acquisition controller 44 that receives signals from a number of 
different sensors connected to the patient, such as ECG signals from electrodes 
attached to the patient. And finally, the pulse generator module 38 connects to a scan 
room interface circuit 46 which receives signals from various sensors associated with the 
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condition of the patient and the magnet system. It is also through the scan room 
interface circuit 46 that a patient positioning system 48 receives commands to move the 
patient to the desired position for the scan. 

[0022] The gradient waveforms produced by the pulse generator module 38 are 
applied to the gradient amplifier system 42 having G x , G y , and G z amplifiers. Each 
gradient amplifier excites a corresponding physical gradient coil in a gradient coil 
assembly generally designated 50 to produce the magnetic field gradients used for 
spatially encoding acquired signals. The gradient coil assembly 50 forms part of a 
magnet assembly 52 which includes a polarizing magnet 54 and a whole-body RF coil 
56. A transceiver module 58 in the system control 32 produces pulses which are 
amplified by an RF amplifier 60 and coupled to the RF coil 56 by a transmit/receive 
switch 62. The resulting signals emitted by the excited nuclei in the patient may be 
sensed by the same RF coil 56 and coupled through the transmit/receive switch 62 to a * 
preamplifier 64. The amplified MR signals are demodulated, filtered, and digitized in the 
receiver section of the transceiver 58. The transmit/receive switch 62 is controlled by a 
signal from the pulse generator module 38 to electrically connect the RF amplifier 60 to 
the coil 56 during the transmit mode and to connect the preamplifier 64 to the coil 56 
during the receive mode. The transmit/receive switch 62 can also enable a separate RF 
coil (for example, a surface coil) to be used in either the transmit or receive mode. 

[0023] The MR signals picked up by the RF coil 56 are digitized by the 
transceiver module 58 and transferred to a memory module 66 in the system control 32. 
A scan is complete when an array of raw k-space data has been acquired in the memory 
module 66. This raw k-space data is rearranged into separate k-space data arrays for 
each image to be reconstructed, and each of these is input to an array processor 68 
which operates to Fourier transform the data into an array of image data. This image 
data is conveyed through the serial link 34 to the computer system 20 where it is stored 
in memory, such as disk storage 28. In response to commands received from the 
operator console 12, this image data may be archived in long term storage, such as on 
the tape drive 30, or it may be further processed by the image processor 22 and 
conveyed to the operator console 12 and presented on the display 16. 
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[0024] The present invention is directed to a method and system of quantitatively 
assessing diffusion in a region-of-interest or tissue. In particular, the present invention is 
well-suited for evaluating diffusion weighted imaging data acquired from a brain to 
identify and delineate between normal and pathological tissue. As will be discussed in 
greater detail, the present invention is directed to a technique of quantitative evaluation 
of diffusion weighted MR imaging data such that an index or metric of diffusion is 
determined and appropriately manipulated to assess diffusion of a substance, such as 
water or a contrast agent, in a tissue. Generally, the metric is defined as: 

f^a n -a){b n -b)\ 
G(A,B) = " =1 ... 1; 



n=l n=\ 



where {a n } and {b n } are diffusion values measured in N different directions for two , 
independent voxels A and B, and a and b are respective average diffusion values for 
the two voxels. This normalized covariance is slightly different than the standard 
definition of correlation, in. which the denominator is the product of the standard 
deviations of {a n } and {b n } rather than their energies. This measure does not assume 
that diffusion is fully characterized by a tensor (i.e. no assumption of a single diffusion 
compartment), and may be used with arbitrary N. Eqn. 1 yields the following 
observations. 

[0025] First, if either voxel A or B is isotropic, the covariance term in the 
numerator (and hence G) will be 0. Second, this covariance term will be positive if the 
larger diffusion terms in A and B are aligned, and smaller (or negative) if they are not. 
Third, the normalization in the denominator keeps G from varying with the average 
diffusion in A or B. Thus G reflects only the shape and relative orientation of the 
measured diffusion in these two voxels. Despite the independence of this metric from 
tensor assumptions, the above observations may be further supported by analytically 
calculating G using six diffusion directions, N=6, with a conventional tensor model. 
According to the prior art tensor-based model, in voxel A, diffusion is oriented along the 
axes {X':Y':Z} (D xr = Dxz = D Y z = 0), where {X':Y'} is rotated from {X:Y} about the Z axis 
by an angle a. Diffusion along X', Y\ and Z is equal to f ax D A) f ay D A , and faz D Al 
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respectively, where D A represents an arbitrary diffusion value and f ax ,y,z represent 
directional deviations from this arbitrary diffusion value. For each direction n, the 
measured diffusion a n for voxel A is equal to: 

a n =D A (f„ cos 2 [a+6IJsin 2 [0J + f ay sin 2 [a +0Jsin 2 [0J + f az cos 2 ty„]), Eqn. 2; 

where O n and 9 n are the spherical angles corresponding to the selected direction relative 
to the +Z axis and the +X axis, respectively. A similar equation is used for voxel B, but 
with diffusion values f bx D B , f by D B , and f bz D B along axes {X":Y M :Z} (rotated by angle (3 
about the Z axis). 



[0026] Inserting the values of a n and b n (n=1 to 6) from Eqn. 2 into Eqn. 1 and 
simplifying yields the following: 

(f, +f ay -2fjf bx +f by - 2f J+3(f. -fjfc, -f by )Cos[2(a-f})] ^ g 

' 3^2 (£ +fi +£) + (f 0 + f ay +f J fiti + f£ + f£)+(f ta +f by + fj 

The angular dependence in Eqn. 3 is only on (a-P), showing that G(A,B) is rotationally 
invariant along the Z axis. It should be noted that the global diffusion values D A and D B 
have dropped out due to normalization. This equation, which is consistent with tensor- 
based modules, is still quite large to be intuitive, so further assumptions may be made. 
For example, assuming voxels A and B have altered diffusion only along the X' and X" 
axes, respectively; the following assignment may be made: f a2 = f ay = fbz = fby = 1 , fax = f a> 
and f bx = f b in Eqn. 3. F 3)b > 1, is consistent with a "cigar" diffusion model and f a , b < 1 is 
consistent with a "pancake" diffusion model. This now yields: 

G(A,B)= , , , 1 + 3C f a) Eqn. 4. 

V(3£+^+8)V(3£+44+8) 3 

In Eqn. 4, a now represents the respective angle between the rotated x-axes of voxels A 

and B. The first two terms in Eqn. 4 are anisotropy terms; they are 0 for f a = 1 or f b = 1 , 

respectively, and monotonically increase with f a b . The last term is an alignment term; it 

is maximum for a = 0° or 180°, and minimum for a = +/- 90°. This analysis illustrates that 

G is (1) large only if the diffusion values in both A and B are similarly anisotropic (e.g. f a 

> 1 and f b > 1 , or f a < 1 and f b < 1), (2) zero if the diffusion in either A or B is isotropic, 
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and (3) largest when the diffusion anisotropics of A and B are 'aligned' (i.e. a = 0° or 
180°). 

[0027] As will be described below, the index or metric defined by Eqn. 1 may be 
manipulated to derivatively assess fractional anisotropy characteristics in a particular 
voxel, filter diffusion data to smooth noise between neighboring voxels, and color-code 
diffusion data. The following are examples illustrating applications in which the present 
invention may be implemented and is not exhaustive of those applications to which the 
invention may be applied as other applications are contemplated and considered within 
the scope of this invention. 

[0028] As defined by Eqn. 1, the value G may be used to rapidly obtain a 
measure of fractional anisotropy with arbitrary N. In particular, the anisotropy metric Gaa 
involves setting all b n = a n in Eqn. 1, and taking the square root of the result, i.e.: 



The numerator in Eqn. 5 becomes the standard deviation of {a n }, and the denominator 
normalizes the metric by the root mean squared diffusion value. The metric in Eqn. 5 
uses the normalized standard deviation of diffusion rather than its variance (i.e. 
A /G(A,A) rather than G(A,A)) so that its behavior is more similar to a standard fractional 
anisotropy measure. 

[0029] In the following example, a comparison between the present invention 
and a prior art, tensor-based approach is described. For N=6, an analytical relationship 
between Gaa and the standard fractional anisotropy measure F can be established as 
follows. Using six selected directions of arbitrary orientation, the tensor parameters may 
be shown to be: 




Eqn. 5. 
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Eqn. 6; 



where c, = 



105 +47 
160 +72V5' 
1 



c 2 = -T. 



c 3 = 



4 

15 +7V5 
160 +72 V5 ' 

45 +20^5 



Eqn. 7; 



4 ~ 80 +3675 ' 

the D vector represents diffusion terms corresponding to the gradient axes, and the a 
vector contains the diffusion values along the directions of the applied diffusion 
gradients. Defining 
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Eqn. 8; 
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Eqn. 9; and 



a = 3 a 2 

mag 
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Eqn. 10. 



The conventional fractional anisotropy metric may be calculated using: 

Eqn. 11. 



■t 



mag 



surface 



Using Eqns. 6-11, an expression for F in terms of a„ may be obtained. While this 
expression is quite long to be of use, this expression may be used along with Eqn. 5 to 
show that: 

2F 



GAA = ViT 



6F 2 



Eqn. 12. 
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One skilled in the art would appreciate that for 0<F<1, this equation is nearly linear. A 
plot of a two-dimensional histogram of F and G M from the same image (including noise), 
would reveal that F and Gaa follow this curve exactly, with no residue. As such, F may 
be obtained directly from Gaa- Moreover, the commonality in the plots would show that 
Gaa and F give equivalent information when N=6 for the selected directions of diffusion. 
Thus, Eqn. 5 illustrates determining of fractional anisotropy for N=6 without the 
determining of a diffusion tensor. Furthermore, Gaa may be calculated for arbitrary N. 
Numerous comparisons of Gaa and F from in-vivo data have confirmed the unique one- 
to-one relationship in Eqn. 5 for every voxel, including noise, when N=6. 

[0030] The present invention is also applicable to nonlinear filtering of DTI data 
prior to tractography, orientational display, or fractional anisotropy calculations. In this 
regard, nonlinearly averaging of every voxel in a 3D data set with itself and M neighbors 
may be carried out using the following expression: 



where p n (B) is the MRI signal value for voxel B for the n th diffusion measurement, p n is 
the new (filtered) data, G(A,B) = max[G(A,B),0], and both summations are for voxel B 
equal to voxel A as well as its M neighbors. Using G(A,B) simply turns all negative 
coefficients to 0; G(A,B) will be negative for crossing, perpendicular tracts, thereby 
allowing for smoothing of data within tracts, rather than let separate tracks affect each 
other's data. 

[0031] From Eqn. 4, it can be expected that the output data p n are significantly 
averaged with that of neighboring voxels if and only if: (a) both voxels are highly (and 
similarly) anisotropic, and (b) the diffusion directions of both voxels are "aligned". Such 
a filter may prevent white matter at the edge of a tract from being smoothed with 
surrounding tissue with isotropic diffusion. The filter also prevents tracts which cross 
each other from being averaged, due to the dependence of G(A,B) on a. The filter may 
decrease anisotropy which may be due to noise reduction, which itself increases 
anisotropy. 




Eqn. 13; 
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[0032] The present invention is also applicable to show directional information in 
a single image. A number of techniques have been developed to generate such a 
display. Most of these techniques are predicated on an "absolute value" method for 
color display, in which the red, green, and blue channels correspond to Dxx, Dyy, and 
Dzz multiplied by an anisotropy index. This can be rapidly approximated according to the 
present invention in the following way. First, N artificial diffusion values for unidirectional 
diffusion along each of the 3 axes (x, y, and z) are synthesized. These values may be 
referred to as {x n }, {y n }, and {z n }. For example, {x n } are synthesized using Eqn. 2 with a = 
0, fxx D x = 1 , and fxy = fxz = 0, leaving the following: 

x n - cos 2 [9J sin 2 [(|) n ] Eqn. 14; 

where 8 n and O n are the spherical coordinates of the N measurement directions. These 
measurements represent unit diffusion oriented purely along the x-axis. 

[0033] Next, for each voxel A, three values (P, T, and B) are calculated as in the 
following example: 

?(A) = G(A,X) +^ Eqn. 15; 

where X represents a "virtual" voxel with diffusion values {x n }. The values r and B are 
calculated (e.g.) using G(A,Y) and G(A,Z) instead of G(A,X), respectively. The value P 
is assigned to the (e.g.) red channel for voxel A, while (e.g.) !~ (A) and B(A) are assigned 
to the green and blue channels, respectively. The rationale for this assignment is that 
the values of P, T, and B are proportional to the anisotropy of A and its orientation to 
each of the 3 axes, respectively. 

[0034] The calculation of P(A) is now illustrated in more detail for N=6, using a 
simplified, conventional tensor model. For voxel A, it is assumed that f ay = f a2 = 1 and f ax 
= f a (i.e. following the assumptions leading up to Eqn. 4 for which the diffusion is 
anisotropic along a single direction x). The above-mentioned values B = 0, D x = 1, fxy = 
f^ = 0 and fxx = 1 may be substituted for voxel B (now the synthesized diffusion along X) 
in Eqn. 4 to get the following expression: 



14 



Docket No. 133625MR(GE 166) 



_ t-i^ft-QGupa] Eq , 16; 

3 V3 V3f a 2 +4f a +8 

^ = - Eqn. 17; 

V3V3f a 2 +4f a +8 



and 



?(A)= 2(f-l)Cos 2 [a] 

V3 V3f a 2 +4f a -f8 Eqn. 18; 

= G AA Cos 2 [aJ 

Thus, with a tensor model and the assumption f ay = f az = 1, P(A) is proportional to the 
anisotropy of A multiplied by a smoothly varying angular function which is 1 when A is 
aligned with the x axis, and 0 when it is perpendicular. Eqn. 18 leads to the relationship 

P(A) + T(A) + B(A) = Gaa Eqn. 1 9, 

i.e. the red, green, and blue channels for any voxel add up to the anisotropy measure of 
the present invention, while the color indicates tract orientation. As such, the present 
invention does not require calculation of diffusion tensor values, nor of eigenvectors. 
The present invention can also be used with an arbitrary number N of diffusion 
measurements, i.e. more than six. 



[0035] The present invention provides is a robust, functional, and rapid method 
for assessing diffusion anisotropy in a single voxel, or shared anisotropy between 
voxels. The techniques of the present invention provide a tensor-free selection that can 
be carried out with an arbitrary number of diffusion directions. The present invention, as 
will be described with respect to Figs. 2-3, includes a method of processing and 
evaluating diffusion weighted imaging data that may be carried out, automatically, as a 
computer program by a one more computers, integrated or independent from an MRI 
system. 

[0036] Referring now to Fig. 2, the steps of a process 70 for determining 
diffusion similarity between two, preferably adjacent, voxels are set forth. The process 
70 begins at 72 with the acquisition of diffusion weighted MR data in accordance with 
one of a number of known data acquisition and imaging techniques. After the MR data 
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is acquired, two adjacent voxels are selected from a region-of-interest, i.e. targeted 
tissue. An average of diffusion throughout each voxel is then determined at 76. Each 
average is determined in accordance with well-known assessment techniques such that 
a value of diffusion is calculated along each direction of diffusion in each voxel. 
Thereafter, individual diffusion values for one or more selected directions are determined 
at 78. The difference between each individual diffusion value and the average difference 
value is then summed at 80. From the individual diffusion values and the average 
diffusion value, a diffusion index or metric G(A,B) as defined by Eqn. 1 is determined at 
82. After G(A,B) is determined for the selected voxels, process 70 either prompts the 
user or automatically determines whether more voxels are to be inspected at 84. If so 
84, 86, process 70 returns to step 74. If not 84,88, the process ends at 90. 

[0037] Regarding G(A,B), if either voxel is isotropic, the covariance term in the 
numerator, as well as the index, will be 0. Also, this covariance term will be positive iP 
the larger diffusion terms in the voxels are aligned, and smaller (or negative) if they are 
not. Further, the normalization in the denominator keeps the index from varying with the 
average diffusion in either voxel. Thus the value of the index will reflect only the shape 
and relative orientation of the measured diffusion in the selected two voxels. 

[0038] Referring now to Fig. 3, the steps of a technique in accordance with 
another embodiment of the invention are step forth. Process 94 is implemented to 
determine fractional anisotropic behavior in a single selected voxel. Process 94 is 
similar the technique of Fig. 2, but is used to quantify diffusion in a single voxel rather 
than multiple adjacent voxels. Process 94 begins at step 96 with the acquisition or 
accessing of diffusion weighted MR data. The data may be newly acquired MR data or 
accessed from a memory location. At step 98, a voxel is selected for evaluation. 
Accordingly, an average diffusion value for the selected voxel is determined or otherwise 
calculated at 100 as is a value of diffusion relative to each selected direction. At step 
102, the process carries out the necessary processing and calculations to solve Eqn. 5 
discussed above. From the output value of Eqn. 5, as described previously, fractional 
anisotropic behavior in the selected voxel may be determined. Process 94 then ends at 
104 or may be repeated to assess fractional anisotropic behavior in another voxel. 
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[0039] The present invention may also be used to develop a filter for smoothing 
noise between adjacent voxels. That is, for two adjacent voxels, A and B, G(A,B), as 
defined by Eqn. 1, is large, if and only if, voxels A and B are anisotropic and "pointing" in 
the same direction. Specifically, G(A,B) may be expressed as: 

G(A,B) = Gaa * G BB * Cos (26) (Eqn. 20); 

where Gaa and G B b correspond to fractional anisotropy wherein voxels A and B, 
respectively, and 9 corresponds to the angular difference between the primary diffusion 
directions in voxels A and B. Further, because of the relationship defined by Eqn. 20, 
the values from Eqn. 1 and Eqn. 5 may be used to determine the relative difference in 
directions. As a result, G(A,B), as determined by Eqn. 1, may be used as a coefficient 
for smoothing neighboring voxels to reduce noise. For example, for brain imaging, only u 
white matter voxels that are in the same tract are averaged together, which preserves 
small tracts and tract edges. In this regard, this embodiment of the present invention is 
particularly well-suited for preparing MR data for tractography processes. 

[0040] Therefore, in accordance with one embodiment, the present invention 
includes a method of assessing diffusion of a substance in a region-of-interest. The 
method includes the steps of determining a diffusion value for one or more individual 
directions in at least one voxel and determining a value of average diffusion in the at 
least one voxel. The method further includes the step of, from the individual diffusion 
values and the average diffusion value, determining a diffusion index indicative of 
diffusion along the first direction in the at least one voxel. 

[0041] According to another embodiment, the present invention includes an MRI 
apparatus comprising a magnetic resonance imaging (MRI) system having a plurality of 
gradient coils positioned about a bore of a magnet to impress a polarizing magnetic field 
and an RF transceiver system and an RF switch controlled by a pulse module to transmit 
RF signals to an RF coil assembly to acquire MR images. The system also includes a 
computer programmed to access diffusion-weighted imaging data from a region-of- 
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interest and isolate at least two adjacent voxels within the region-of-interest. The 
computer is further programmed to determine a metric indicative of diffusion along 
multiple directions in the at least two voxels without fitting diffusion-weighted imaging 
data from the at least two adjacent voxels to a tensor model. 

[0042] In accordance with another embodiment of the present invention, a 
computer readable storage medium having stored thereon a computer program 
representing a set of instructions is provided. The set of instructions, when executed by 
a computer, causes the computer to determine an individual diffusion value for each of N 
diffusion directions in at least one voxel and determine an average diffusion value for the 
at least one voxel. The set of instructions further causes the computer to determine a 
diffusion index indicative of diffusion anisotropy and diffusion direction in the at least one 
voxel. 

[0043] The present invention has been described in terms of the preferred 
embodiment, and it is recognized that equivalents, alternatives, and modifications, aside 
from those expressly stated, are possible and within the scope of the appending claims. 
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